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INTRODUCTION 


When  a  high-velocity  projectile  impacts  and  penetrates  a  semi-infinite  body  of  con¬ 
crete,  the  penetration  is  accomplished  by  displacing  the  target  material  radially  and  as  a 
result  a  tunnel-shaped  crater  is  formed  in  the  target.  The  target  resistance  to  penetration 
(which  is  inferred  by  the  depth  and  diameter  of  the  produced  hole)  depends  on  a  number  of 
parameters  including  the  projectile  velocity  and  the  relative  strength  of  the  projectile  and 
the  target.  In  particular,  when  the  penetrator  velocity  is  fixed,  the  structure  of  the  plastic 
flow  field  in  the  target  defines  the  extent  of  the  ensuing  target  material  displacements  (both 
the  radial  and  the  axial)  which  in  turn  determines  the  size  of  the  resulting  crater.  Therefore, 
since  the  structure  of  this  flow  field  is  controlled  primarily  by  the  yield-strength  properties  of 
the  target,  the  constitutive  behavior  of  the  target  material  is  the  principle  factor  that  deter¬ 
mines  the  target's  resistance  to  penetration. 

Under  the  pressures  exerted  at  the  interface  between  the  projectile  and  the  target,  the 
penetrator  nose  may  deform  into  a  characteristic  mushroom-like  shape.  The  extent  of  this 
deformation  varies  significantly  with  the  impact  velocity  and  the  relative  strength  of  the  pro¬ 
jectile  and  the  target.  For  impact  velocities  exceeding  1 .2  to  1 .5  km/s,  the  strength  of  most 
penetrator  materials  is  not  sufficient  to  prevent  penetrator  deformation  and  mushrooming, 
resulting  in  significant  penetrator  erosion.  Recent  experimental  work  by  Gold  et  al.  (ref  1 ) 
showed  that  for  velocities  exceeding  1 .5  km/s  the  penetrator  erodes  almost  entirely,  im¬ 
plying  that  penetration  of  a  semi-infinite  body  of  concrete  can  be  successfully  predicted 
employing  the  one-dimensional  modified  (i.e.,  strength  dependent)  hydrodynamic  theory  of 
penetration  (ref  2).  Within  the  framework  of  this  theory,  the  essence  of  the  mechanics  of 
penetration  of  a  solid  projectile  into  a  solid  target  is  represented  through  erosion  of  the 
penetrator  material  in  combination  with  its  rigid  body  motion.  According  to  the  theory,  the 
strength  properties  of  the  penetrator  and  the  target  are  represented  with  two  empirical 
strength  factors  Ypand  Rt,  usually  referred  to  as  the  penetrator  and  target  "strengths," 
respectively.  For  each  combination  of  the  penetrator  and  target  materials,  the  empirical 
strength  factors  Ypand  Rt  have  to  be  determined  experimentally.  One  of  the  principle  con¬ 
clusions  of  the  modified  hydrodynamic  theory  of  penetration  is  that  the  target  resistance  to 
penetration  is  proportional  to  the  target  strength  factor  Rt,  a  fact  supported  by  a  vast 
amount  of  experimental  data. 

The  analytical  structure  of  the  parameter  Rt  has  been  examined  by  a  number  of 
researches  including  Tate  (refs  3  and  4),  Pidsley  (ref  5),  and  Anderson  and  Walker  (refs  6 
and  7)  both  analytically  and  numerically.  Because  the  assumptions  of  the  modified  hydro- 
dynamic  theory  of  penetration  neglect  the  complexity  of  the  structure  of  the  target  elastic- 
plastic  flow  field,  the  discrepancy  between  the  experimental  values  of  Rtand  the  numerical 
predictions  is  rather  significant.  For  example,  in  contradiction  to  the  assumptions  of  the 
modified  hydrodynamic  theory,  the  numerical  results  of  Anderson  and  Walker  (refs  6  and  7) 
showed  that  during  the  penetration  the  magnitude  of  Rt-Vp  is  not  constant,  but  changes  as 
the  penetrator  penetrates  the  target.  The  analysis  presented  in  this  work  examines  the 
relationship  between  the  target  strength  factor  Rt  and  the  target  resistance  to  penetration. 
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Concrete  is  a  complex  material  consisting  of  mineral  aggregate  bound  by  cement 
paste  containing  a  large  amount  of  water  and  voids.  In  the  present  work,  two  concrete 
yield-strength  models  were  studied:  (1)  the  constant  yield-strength  model  and  (2)  the 
pressure-dependent  yield-strength  model.  In  the  case  of  the  constant  yield  model,  the 
disagreement  between  the  calculated  hole  profiles  and  the  experimental  data  was  rather 
significant,  regardless  of  the  value  of  the  yield  strength  parameter  k.  In  the  case  of  the 
pressure-dependent  yield-strength  model,  the  analysis  resulted  in  hole  predictions  which 
not  only  reproduced  the  experimental  depth  of  penetration  but  also  matched  the  measure¬ 
ments  of  the  profile  of  the  tunnel  portion  of  the  crater. 

The  reduction  in  the  calculated  penetration  channel  diameters  (with  the  introduction  of 
the  pressure-dependent  yield  strength)  was  attributed  to  the  increased  target  resistance  to 
penetration.  In  order  to  clarify  the  effect  of  the  constitutive  modeling  on  the  penetration  cal¬ 
culations,  further  analysis  concentrated  on  studying  the  parameters  of  the  flow  field  along 
the  central  streamline.  Along  the  centerline  a  number  of  flow  field  parameters  vanish,  thus 
permitting  significant  simplification  of  the  momentum  balance  equation.  Assuming  that  the 
penetration  is  steady  state,  the  momentum  balance  equation  can  be  integrated  resulting  in 
an  explicit  expression  for  the  pressure  at  the  target/penetrator  interface.  In  this  expression, 
the  target  resistance  to  penetration  is  represented  by  the  following  terms:  (1)  the  difference 
between  stress  deviators  at  the  elastic-plastic  boundary  and  the  target/penetrator  interface, 
respectively,  (2)  the  pressure  at  the  elastic-plastic  boundary,  and  (3)  the  integral  of  the 
radial  gradient  of  the  shear  stress  deviator  (at  the  central  streamline)  evaluated  within  the 
limits  of  the  plastic  zone.  A  detailed  analysis  of  the  above  terms  shows  that  inclusion  of  the 
yield  strength-pressure  dependency  in  the  constitutive  model  significantly  increases  the 
target  resistance  to  penetration.  As  a  result,  the  time  rate  of  projectile  erosion  increases, 
which  ultimately  leads  to  significant  reduction  in  the  calculated  penetration  depths  and  hole 
diameters. 


CONCRETE  PENETRATION  EXPERIMENTS  AND  THE  NUMERICAL  MODEL 

The  experiments  and  numerical  modeling  referred  to  in  the  present  study  have  been 
reported  extensively  in  the  literature  (refs  1  and  8)  and  will  be  briefly  presented  here  for 
completion.  Copper  and  tantalum  projectiles  were  launched  by  gun  against  semi-infinite 
concrete  and  reinforced  concrete  targets  with  velocities  ranging  from  0.15  cm/ps  to  0.19 
cm/ps.  The  projectiles  were  spherical-nose,  right  circular  cylinders  with  diameters  of 
1 .3  cm  and  2.0  cm,  and  length-to-diameter  ratios  varying  from  3.9  to  14.6.  The  projectile's 
launching  velocities  were  measured  from  streak  camera  records.  Orthogonal  flash  x-rays 
in  the  immediate  vicinity  of  the  surface  of  the  target  independently  verified  the  impact 
velocities  and  provided  information  on  the  yaw  and  condition  of  projectiles  just  prior  to  the 
impact. 
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The  targets  were  right  circular  cylinders  91  cm  in  diameter  and  91  cm  long  and  were 
constructed  from  concrete  with  maximum  aggregate  size  of  1 .9  cm  and  density  of  2.24 
g/cm^.  Laboratory  testing  of  specimens  of  this  concrete  resulted  in  an  average  unconfined 
compressive  strength  of  approximately  0.374  Kbar.  The  targets  were  instrumented  with  a 
sequence  of  break  gages  that  were  carefully  positioned  with  reference  to  the  front  surface 
and  spaced  2.54  cm  apart  from  each  other.  The  gages  were  0.03-cm  thick  printed  circuit 
boards  consisting  of  two  conductive  copper  foils  and  enclosing  an  insulated  copper  foil 
maze  between  them.  Although  it  is  rather  difficult  to  define  the  exact  mode  of  gage  failure, 
it  is  reasonable  to  assume  that  once  the  strain  exceeds  that  of  the  elastic  limit  of  the  weak¬ 
est  material  (insulation),  the  gage  "breaks."  When  the  projectile  penetration  rate  is  below 
the  sonic  velocity  of  the  target  material,  concrete  that  is  crushed  in  front  of  the  projectiles 
forms  a  well  defined  plastic  zone  of  a  finite  size.  Therefore,  assuming  that  the  "strengths" 
of  the  gage  and  the  concrete  are  approximately  the  same,  the  gage  would  "break"  approxi¬ 
mately  at  the  time  when  the  elastic-plastic  boundary  arrives  at  the  position  of  the  gage. 
Thus,  the  succession  of  trigger  times  of  the  break  gages  determines  the  trajectory  of  the 
elastic-plastic  boundary  in  front  of  the  projectile. 

In  all  experiments  the  resulting  crater  profiles  were  carefully  measured.  In  the  case  of 
the  1.3-cm  diameter  projectiles,  crater  entrances  were  approximately  25  cm  in  diameter;  the 
penetration  channel  rapidly  narrowed  down  to  a  diameter  of  5  cm,  10.5  cm  away  from  the 
surface  of  the  target.  The  deeper  portion  of  the  crater  had  a  characteristic  form  of  a  well 
rounded  and  slightly  tapered  tunnel.  Measurements  of  this  portion  of  the  crater  were  used 
for  validating  the  accuracy  of  different  constitutive  models  for  concrete  studied  in  this  work. 

In  the  present  work,  the  numerical  experiments  concentrated  on  modeling  a  repre¬ 
sentative  experiment  with  a  19-cm  long  and  1.3-cm  diameter  copper  projectile  impacting  a 
plain  concrete  target  with  velocity  \/o=0.1836  cm/ps  [test  PA8Cu1.4  (ref  1)].  The  numerical 
analysis  presented  here  was  performed  using  the  CALE  (ref  9)  computer  program.  CALE  is 
a  plane  two-dimensional  and  three-dimensional  axisymmetric  hydrodynamics  code  based 
on  an  Arbitrary  Lagrangian  Eulerian  formulation  of  the  governing  equations.  It  can  run  in  a 
pure  Lagrangian  mode,  a  pure  Eulerian  mode,  or  in  an  arbitrary  hybrid  mode.  In  the 
Eulerian  reference  frame  the  description  of  the  axisymmetric  kinematic  deformation  of  the 
continuum  is  given  by  the  following  conservation  laws: 

Conservation  of  mass 

dp/dt  +  Vrdp/dr+  Vrdp/dz  +  p{{Mr)d{rv,)ldr  +  dVj/dz)  =  0  (1 ) 

Conservation  of  linear  momentum 

dVrIdt  +  VrdVf/dr  +  Vz-dvjdz  =  {;Mp){doJdr  +  d{on^p)/dz  +  {on+  Ozz+2p)/r))  (2a) 

dvjdt  +  VrdVzfdr+  Vrdvjdz  =  (1  /p)(aaa/az  +  a(a„+p)/ar  +  (oK+p)/r)  (2b) 
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Conservation  of  energy 

dE/dt  +  v,  dE/dr-^  Vz  dE/dz  =  {Mp)  [-p  {(Mr)  d{r  Vrydr  +  dVg/dz)  +  r  {d(Vrlr)ldr)  (a„+p) 

+  {aa+p){dvzldz-vjr)  +  {cin^-p){dVjldz-dvJdr)\  (3) 

Bcxly  forces  (such  as  gravity)  and  heat  conduction  are  ignored. 

GALE'S  procedure  for  calculating  elastic-plastic  flow  is  based  on  the  methodology 
developed  by  Wilkins  (refs  10  and  1 1 ).  The  stress  behavior  of  a  material  can  be  thought  of 
being  composed  of  a  stress  associated  with  uniform  hydrostatic  pressure  p,  referred  to  as 
the  spherical  part  of  the  stress  tensor,  plus  the  stress  associated  with  the  resistance  of  the 
material  to  shear  distortion,  referred  to  as  the  deviator  Sy.  In  describing  yielding  and  plastic 
flow,  the  stress  contributions  that  are  due  to  shear  distortion  are  limited  by  the  constitutive 
equation  for  the  material.  The  stress  tensor  is  given  then  by 

Oy=-p  +  Sy  (4) 

The  constitutive  equation  for  the  projectile  material  adopted  here  is  that  of  Steinberg- 
Guinan  (ref  12),  which  specifies  the  shear  modulus  and  the  yield  strength  as  a  function  of 
pressure,  temperature,  and  equivalent  plastic  strain.  The  equation  of  state  (EOS)  that  was 
used  to  model  the  hydrodynamic  response  of  copper  was  the  standard  linear  polynomial 
approximation  usually  employed  for  metals  and  adopted  by  Tipton  (ref  9)  for  the  GALE 
code.  Since  the  copper  projectiles  were  manufactured  from  a  standard  stock  of  oxygen- 
free  high  conductivity  (OFHG)  copper,  the  analysis  relied  on  a  standard  set  of  parameters 
for  this  material  given  by  Tipton  (ref  1 3).  Details  of  these  models  and  their  implementation 
can  be  found  in  Tipton  (ref  9). 

In  modeling  the  constitutive  response  of  concrete  targets,  all  models  studied  account¬ 
ed  for  compressibility  of  conaete.  The  nonlinearity  in  the  Hugoniot  data  for  the  concrete 
for  particle  velocities  below  1  km/s  has  been  attributed  to  the  effects  of  microcracking  and 
pore  closure,  and  was  modeled  with  a  porous  material  equation  of  state.  Using  a  porous 
EOS  model.  Read  and  Maiden  (ref  14)  were  able  to  reproduce  the  essential  features  of  the 
experimental  shock-velocity  versus  particle-(material)-velocity  data.  The  particle  velocities 
experienced  by  the  concrete  in  the  earlier  penetration  experiments  of  Gold  et  al.  (ref  1 )  are 
in  the  very  same  range.  Therefore,  employment  of  the  porous  EOS  model  is  ideally  suited 
to  the  conditions  of  high  velocity  penetration. 

Following  Tipton  (ref  15),  the  porous  EOS  model  assumes  that,  while  in  compression, 
the  porosity  is  not  changed  until  the  pressure  exceeds  the  compaction  pressure,  Pc.  In 
tension,  the  porosity  is  not  increased  until  the  pressure  drops  below  -pc.  The  compaction 
pressure  is  assumed  to  be  a  function  of  the  porosity  cp  and  has  the  following  form 


Pc=-Po/og((p) 


(5) 
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The  total  pressure  of  the  porous  material  is  given  by 


p  =  (1-(p)Ps(Ps)  (6) 

where  cp  is  the  porosity  and  psand  psare  the  pressure  and  the  density  of  the  solid  phase, 
respectively.  The  EOS  for  the  solid  phase  is  taken  in  a  polynomial  form 

Ps=Ao  +  /Aip  +  A2[i^  +  +  (So  +  Bip  +  B2[i^)E  (7) 

where  E  is  the  internal  energy,  p=p/pr1 ,  pr=po/(1-^o)  is  the  reference  density,  and  (po  is  the 
initial  porosity.  Further  details  on  this  model  and  its  implementation  in  GALE  can  be  found 
in  Tipton  (ref  1 ). 

GALE'S  procedure  for  calculation  of  the  elastic-plastic  flow  employs  the  von  Mises 
yielding  criterion.  The  yielding  behavior  of  concrete  was  modeled  using  the  following 
constitutive  equations: 


(1)  Constant  yield-strength  model 

J2=  l^=  Const 

(2)  Pressure  dependent  yield-strength  model 

J2=f^{P) 

In  these  models,  J2  is  the  second  deviatoric  stress  invariant  given  by 

>-^2  ~  ~  ■  (Srr  Sq,(p  +  S^z  )  Sfcp  "*■  Srz  ■*"  Szjp 


(8) 

(9) 


(10) 


and  k  or  k{p)  are  either  an  empirical  strength  constant  or  an  empirical  function  (of  pressure 
p),  respectively.  The  relationship  between  the  experimental  strength  data  and  k  is  as 
follows:  consider  an  ideal  test  in  which  a  specimen  is  subject  to  a  simple  tension/compres¬ 
sion  in  the  z-direction  (i.e.  Ozz=Y),  while  the  rest  of  the  stress  components  vanish  Oij=0.  In 
the  case  of  a  constant  strength  model  given  by  equation  8,  J2=Y^I3,  or  k=YIY^.  In  the  case 
of  the  pressure  dependent  yield  model,  the  empirical  strength  function  k{p)  is  obtained 
directly  from  triaxial  compression  test  data  (ref  16). 

The  geometry  and  boundary  conditions  are  shown  in  figure  1 .  The  domain  consider¬ 
ed  is  axisymmetric  in  r-z  space.  Domain  7i  727374  represents  the  target,  P1P2P371  is  the 
spherical-nose  projectile,  and  the  rest  of  the  domain  EiE2727i  is  filled  with  air.  Eight  0.03- 
cm  thick  break  gages  were  initially  modeled  using  three  computational  cells,  each  0.01 -cm 
thick,  filled  with  copper.  After  some  initial  numerical  experimentation  it  became  clear  that 


5 


these  0.01  cm  cells  were  resulting  in  very  small  time  steps,  which  controlled  the  overall 
speed  of  the  calculations.  Therefore,  in  order  to  decrease  the  otherwise  prohibitively 
excessive  computational  time,  the  number  of  cells  across  the  thickness  of  the  break  gages 
was  reduced  to  one.  A  more  detailed  description  of  the  computational  mesh  employed  in 
the  model  can  be  found  in  Gold  et  al.  (ref  8). 

The  initial  and  the  boundary  conditions  are  as  follows: 
at  f=0  Vr=0,  Vz=Vo  in  domain  P1P2P3T1 

Vr=V2=0  everywhere  inside  of  £1527374,  except  P^P2P^T^  (1 1 ) 


at  t>0 

11 

11 

0 

along  £i52 

an=0 

along  £273  and  T3T4 

Vr=0 

along  £1 T4  (line  of  axial  symmetry) 

(12) 


RESULTS  OF  THE  NUMERICAL  ANALYSIS 

The  analysis  presented  in  this  work  addresses  a  number  of  issues  pertaining  to  con¬ 
stitutive  modeling  of  concrete  in  high  velocity  penetration  analysis.  Two  concrete  yield- 
strength  models  were  studied:  (1)  the  constant  yield-strength  model  given  by  equation  8, 
and  (2)  the  pressure  dependent  yield-strength  model,  given  by  equation  9.  Hole  profiles 
resulting  from  the  numerical  experimentation  with  these  models  are  shown  in  figure  2a  and 
b.  In  the  case  of  the  constant  yield  model,  the  disagreement  between  the  calculated  hole 
profiles  and  the  experimental  data  was  rather  significant  (fig.  2a),  regardless  of  the  value  of 
the  yield  strength  parameter  k  [the  von  Mises  yield  strength  k=yi^  was  varied  in  the  rela¬ 
tively  wide  range  from  y=0.187  Kbar  to  y=1 .8  Kbar,  while  the  unconfined  compressive 
strength  of  the  concrete  was  approximately  0.374  Kbar  (ref  1 )].  Since  the  numerical  experi¬ 
ments  resulted  in  almost  identical  penetration  channel  diameters  (which  all  significantly 
disagreed  with  the  experimental  data),  it  follows  then  that  the  poor  predictive  capability  of 
the  constant  yield-strength  model  cannot  be  improved  by  mere  "adjustment"  of  the  strength 
parameter  k. 

In  the  case  of  the  pressure-dependent  yield-strength  model  (eq.  9)  the  analysis 
resulted  in  hole  predictions,  which  not  only  reproduced  the  experimental  depth  of  penetra¬ 
tion,  but  also  matched  the  measurements  of  the  profile  of  the  tunnel  portion  of  the  crater 
(fig.  2b).  It  follows  then  that  mere  inclusion  of  the  yield-strength-pressure  dependency  in 
the  concrete  constitutive  model  results  in  a  rather  significant  decrease  in  the  hole  profile 
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predictions,  apparently  due  to  the  increased  target  resistance  to  penetration.  Thus,  before 
one  even  attempts  to  understand  the  relationship  between  the  target  constitutive  modeling 
and  the  target  resistance  to  penetration,  one  has  to  identify  principle  factors  which  deter¬ 
mine  the  capability  of  a  target  material  to  resist  penetration. 

It  can  be  hypothesize  that  one  of  the  principle  parameters  that  affect  target  resistance 
to  penetration  is  the  pressure  at  the  target/penetrator  interface.  The  calculated  histories  of 
pressure  at  the  target/penetrator  interface  (at  the  centerline  p=0)  are  shown  in  figure  3. 
According  to  numerous  experimental  and  analytical  studies,  the  penetration  process  can  be 
differentiated  into  the  following  three  stages  (e.g,  refs  17  and  18):  (1)  the  initial  transient 
stage.  (2)  the  quasi  steady-state  penetration,  and  (3)  penetrator  deceleration.  As  seen 
from  figure  3,  during  the  initial  transient  stage  (which  lasts  up  to  approximately  40  ps)  and 
during  the  steady-state  penetration  (that  continues  to  approximately  280  ps),  all  the  calcu¬ 
lated  pressure  histories  are  almost  identical.  It  is  only  in  the  final  stage  of  the  penetration, 
when  the  projectile  debris  start  to  decelerate  rapidly  (which  occurs  after  approximately  280 
ps),  that  the  pressure  history  curves  become  distinguishable  from  each  other.  As  seen 
from  figure  3,  these  later  portions  of  the  curves  are  distinguishable  from  each  other  only 
because  of  the  difference  in  their  respective  starting  times  (which  correlate  well  with  hypo¬ 
thetical  changes  in  target  "strengths").  After  "branching  off",  all  curves  are  approximately 
parallel.  This  indicates  that  changes  in  the  constitutive  model  affect  only  the  commence¬ 
ment  of  the  last  stage  of  the  penetration,  while  the  rates  of  the  deceleration  stay  approxi¬ 
mately  the  same.  Therefore,  if  there  exist  a  difference  in  the  penetration  characteristics, 
this  difference  must  manifest  itself  during  the  stage  of  the  steady-state  penetration. 

During  the  quasi  steady-state  penetration,  the  motion  of  the  penetrator  can  be  viewed 
within  the  framework  of  the  one-dimensional  modified  hydrodynamic  theory  of  penetration 
(refs  2,  3,  and  19).  According  to  this  theory,  pressure  at  the  target/penetrator  interface  is 
given  by  the  following 

j  Ppoiv-uf  +  Vp  =  j  -ptou^  +Rt  (13) 

where  ppo  and  pw  are  densities  of  the  penetrator  and  the  target,  respectively.  The  velocity  of 
the  underformed  portion  of  the  penetrator  is  v,  u  is  the  penetration  rate,  and  Yp  and  Rt  are 
empirical  parameters  relating  to  strengths  of  penetrator  and  target  materials,  respectively. 

In  this  simplistic  one-dimensional  model,  the  pressure  at  the  target/penetrator  inter¬ 
face  is  determined  only  as  a  function  of  local  flow  field  variables  such  as  the  velocity,  the 
physical  properties,  etc.  However,  in  the  case  of  the  three-dimensional  analysis,  the  pres¬ 
sure  is  also  influenced  by  the  presence  of  the  free  boundary,  which  reflects  the  traveling 
stress  waves  back  into  the  target  domain.  Inside  the  target  domain  the  pressure  wavelets 
interfere  and  superimpose,  resulting  in  large  fluctuations  of  the  flow  field  parameters  that 
are  sensitive  to  the  continuum  stress  wave  propagation  properties.  Therefore,  in  order  to 
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compare  these  typically  erratic  and  rapidly  changing  parameters,  one  has  to  redefine  them 
on  a  time  average  basis.  Adopting  an  approach  suggested  by  Anderson  and  Walker  (ref 
6),  the  time  averaged  pressure  can  be  defined  as  follows: 


P  = 


(14) 


where  t  is  the  time,  p=p(f)  is  the  pressure  at  the  target/penetrator  interface,  and  7  is  the 
duration  of  the  quasi  steady-state  stage.  However,  even  during  this  quasi  steady-state 
stage  of  penetration,  the  penetrator  slightly  decelerates,  and,  therefore,  the  pressure  at  the 
target/penetrator  interface  changes  with  time.  Since  equation  14  completely  ignores  these 
effects,  the  notion  of  the  time  averaged  pressure  can  be  modified  assuming  a  linear 
dependency  between  the  average  pressure  p’  and  the  time 

p*(f)  =  Po  +  pif  (15) 

Coefficients  po  and  can  be  easily  determined  by  minimizing  the  following  expression 

N 

x\po,  Pi)  =  Z 


Pi  -Po-P^ii 


(16) 


where  piare  the  values  of  pressure  at  given  times  f,,  /=1  ,...,  N,  and  o/is  the  standard 
deviation.  Employing  this  notion  of  the  time  averaged  pressure,  one  can  readily  examine 
the  effect  of  the  target  constitutive  modeling  on  the  pressure  at  the  target/penetrator 
interface.  The  histories  of  the  time  averaged  pressure  p*  at  the  target/penetrator  interface 
are  given  in  the  enlarged  segment  of  figure  3.  These  results  show  that,  in  agreement  with 
the  hydrodynamic  theory  of  penetration,  the  target/penetrator  interface  pressure  increases 
with  increases  in  the  "target  strength,"  a  tendency  that  may  be  intuitively  associated  with 
the  effect  of  the  von  Mises  yield  strength  k  or  /c(p). 

The  empirical  yield  strength  function  (ref  16)  Y-Y{p),  which  was  employed  in  the 
analysis  using  the  pressure-dependent  yield  model,  is  shown  in  figure  4.  During  the  quasi 
steady-state  penetration  the  average  pressure  at  the  target/penetrator  interface  is 
approximately  0.02  Mbar.  Therefore,  according  to  figure  4,  the  corresponding  averaged 
yield  strength  at  the  projectile/target  interface  is  approximately  1 .6  Kbar.  Intuitively,  one 
would  expect  that  the  higher  values  of  the  "average"  strength  have  to  result  in  the  higher 
values  of  the  pressures.  And  yet,  the  calculations  show  (see  p’  in  fig.  3)  that  in  the  case  of 
the  pressure-dependent  yield  model  (the  "average"  yield  is  approximately  1.6  Kbar),  the 
time-averaged  pressure  p  is  higher  than  the  pressure  resulting  from  the  1 .8  Kbar  constant 
yield  model.  It  follows  then,  that  the  target  resistance  to  penetration  is  not  related  simply 
either  to  the  pressure  or  to  the  yield  strength  alone.  In  order  to  identify  the  principle  factors 
affecting  the  target  resistance  to  penetration,  one  has  to  examine  the  structure  of  the 
plastic  flow  field  in  the  target. 
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The  structure  of  the  elastic-plastic  flow  field  in  the  target  can  be  examined  from  the 
results  presented  in  figure  5.  Isolines  of  the  equivalent  plastic  strain  shown  in  figure  5. 
represent  the  approximate  position  of  the  elastic-plastic  boundary  in  the  target.  The 
equivalent  plastic  strain  is  defined  as 


Y 


(17) 


where  s^p,  S2P,  and  sspare  the  principal  plastic  strain  components.  The  equivalent  plastic 
strain  is  one  of  the  principal  parameters  of  the  plastic  flow  field  which  reflects  the  extent  of 
plastic  state,  and  the  isolines  of  y  characterize  the  structure  of  the  deformation  field.  In  the 
elastic  region  the  plastic  strain  components  vanish,  and  y=0.  In  the  plastic  region  y>0, 
whereas  the  isolines  of  a  small  y  mark  the  approximate  position  of  the  elastic-plastic 
boundary. 

Employing  numerical  experimentation,  one  may  attempt  to  construct  a  reasonable 
approximation  to  the  plastic  flow  around  the  penetrator,  an  approach  useful  in  a  number  of 
the  approximate  theories  of  high  velocity  penetration.  It  is  interesting  to  note  that  all 
numerical  experiments  resulted  in  isolines  with  shapes  that  are  approximately  hemispher¬ 
ical,  regardless  of  the  employed  target  constitutive  model.  The  hemispherical  symmetry  in 
the  target  deformation  field  appears  to  be  an  intrinsic  characteristic  of  the  high  velocity 
penetration  process  proper,  rather  than  a  result  of  a  particular  type  of  a  constitutive  model. 
This  observation  is  also  supported  by  the  results  of  high  velocity  penetration  analysis  for 
various  metallic  targets  (e.g.,  refs  5  and  20  through  23). 

Since  the  target  deformation  field  is  hemispherically  symmetric,  the  center  line  r=0 
seems  to  be  the  most  attractive  choice  for  referencing  both  the  motion  of  the  penetrator  in 
the  target  and  the  flow  of  the  target  material  around  the  penetrator.  Referring  to  figure  5,  Zo 
denotes  the  position  of  the  interface  between  the  penetrator  and  the  target,  while  the 
coordinate  of  the  elastic-plastic  boundary  is  given  by  Zep.  Thus,  the  plastic  zone  in  front  of 
the  penetrator  is  described  by  a  single  parameter  A=Zep-Zo.  The  interrelationship  between  A 
and  the  target  resistance  to  penetration  is  as  follows. 

Consider  a  balance  of  energy  for  the  plastic  zone  Zep(f)^2^o(0  enclosed  within  an 
idealized  control  volume,  which  moves  through  the  target  with  a  velocity  equal  to  the  pene¬ 
tration  rate  u.  In  the  frame  of  reference  which  moves  with  velocity  u,  the  energy  (in  the  form 
of  plastic  work)  is  influxed  into  the  control  volume  through  the  boundary  Zq.  The  amount  of 
this  energy  Eo  is  determined  by  the  penetrator  erosion  rate  v-u,  and  ultimately,  by  the 
kinetic  energy  of  the  penetrator.  The  efflux  of  energy  out  of  the  control  volume  Eep  is 
accomplished  (in  the  form  of  elastic  work  on  the  target  material)  through  the  elastic-plastic 
boundary  Zep.  Assuming  that  the  elastic  work  is  negligible  compared  to  the  kinetic  energy  of 
the  penetrator  (Eep»0),  the  entire  kinetic  energy  of  the  penetrator  is  deposited  into  the 
plastic  work  inside  the  control  volume,  which  eventually  results  in  generation  of  heat. 
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Inside  the  control  volume,  the  distribution  of  the  energy  per  unit  volume  is  a  function  of  the 
structure  of  the  plastic  flow  field.  As  a  first  approximation,  the  total  energy  dissipated  inside 
the  control  volume  can  be  viewed  as  the  volume  of  the  plastic  zone  (~A®)  times  a  coefficient 
K.  The  value  of  this  coefficient  characterizes  the  average  dissipation  of  the  energy  per  unit 
volume,  and  is  a  function  of  the  structure  of  the  plastic  flow  field  and  the  target's  strength. 

It  follows  then  that  since  Eo=k-A^  (and  since  the  energy  Eo  is  fixed  by  the  kinetic  energy  of 
the  penetrator),  the  size  of  the  plastic  zone  A  becomes  one  of  the  principle  parameters  that 
determine  the  resistance  of  the  target  to  penetration. 

The  numerical  experiments  revealed  that,  depending  on  the  particular  structure  of  the 
employed  constitutive  model,  the  size  of  the  plastic  zone  A  was  either  significantly  affected 
by  the  model,  or  the  effect  was  negligible.  In  particular,  all  numerical  experiments  with  the 
constant  yield-strength  model  resulted  in  approximately  the  same  value  of  A«1.3  cm,  re¬ 
gardless  of  the  value  of  the  yield  strength  parameter  Y.  In  contrast,  a  comparison  between 
plots  resulting  from  the  constant  yield-strength  model  (A«1 .3  cm)  and  pressure  dependent 
yield-strength  model  (A«1 .7  cm)  showed  that  a  mere  change  in  the  target  constitutive  model 
resulted  in  a  very  significant  30%  change  in  the  size  of  the  plastic  zone. 

In  the  previous  approach,  the  target  resistance  to  penetration  is  associated  with  the 
amount  of  energy  per  unit  volume  dissipated  in  the  plastic  flow.  If  the  target  is  an  ideal 
fluid,  then,  according  to  the  ideal  hydrodynamic  theory  of  penetration  (refs  24  and  25)  the 
penetration  resistance  of  the  target  is  a  function  of  the  target  density  (i.e.,  inertia)  only.  If 
the  target  is  a  solid  with  a  finite  strength,  the  stronger  the  target  the  larger  there  must  be 
the  value  of  the  energy  dissipation  per  unit  volume.  In  contradiction  to  this  intuitive  con¬ 
clusion,  the  numerical  experiments  indicated  the  opposite:  the  higher  the  target  resistance 
to  penetration  was,  the  larger  was  the  size  of  the  plastic  zone  A,  and  therefore  the  average 
of  the  dissipated  energy  per  unit  volume,  K=Eo/A^  was  lesser.  Explanation  of  these  "contro¬ 
versial"  results  needs  a  detailed  analysis  of  the  notion  of  the  "target  resistance  to  pene¬ 
tration." 


ANALYSIS  OF  THE  AXIAL  MOMENTUM  BALANCE  EQUATION 
ON  THE  LINE  OF  AXIAL  SYMMETRY 


Assuming  that  the  penetration  is  approximately  steady  state,  Pidsley  (ref  5),  and  later 
Anderson  and  Walker  (ref  7),  showed  that  on  the  center  line  r=0  the  axial  momentum 
balance  equation  2b  can  be  reduced  to  the  following  equation 


P  Vz 


—  ~T'  (Szz  ■  P)  + 
dz 


r=0 


r=0,  zo-/ <  z  <  Zep  (18) 


where  /  is  the  length  of  the  penetrator  and  Zo  and  Zep  are  axial  coordinates  of  the  target/ 
penetrator  interface  and  the  elastic-plastic  boundary,  respectively.  Consider  an  arbitrary 
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domain  in  the  target  (z  <Zep),  in  which  p=pt(z).  The  pressure  at  the  target/penetrator 
interface  can  be  found  by  integrating  equation  18  over  the  interval  [zo,z*]: 


P(zo)  =  +  p(z*)  +  s^2(zo)  -  sj,z )  -  2  • 


dz 


r=0 


R. 


(19) 


The  term  Rh  appearing  on  the  right  hand  side  of  equation  19  represents  the  "hydro- 
dynamic"  component  of  the  pressure  at  the  target/penetrator  interface,  and  assuming  that 
the  target  material  is  incompressible  (i.e.,  pt=const),  the  integral  can  be  evaluated  easily. 

In  this  work,  concrete  was  modeled  as  a  porous  or  "crushable"  material,  which  has  a 
significant  "built-in"  degree  of  compressibility.  Based  on  the  results  of  numerical  experi¬ 
ments,  it  can  be  shown  that  the  assumption  of  incompressibility  can  also  be  applied  to  this 
particular  material  model.  In  all  cases  studied,  concrete  was  modeled  employing  the 
porous  equation  of  state  given  by  equations  5  to  7  in  which  the  degree  of  the  compressi¬ 
bility  of  this  material  is  controlled  primarily  by  the  compaction  pressure  criterion,  equation  5. 
Since  the  value  of  the  concrete  compaction  pressure  parameter  Pc  [pc=1 .8  Kbar  (ref  1 5)] 
was  greater  than  the  value  of  the  von  Mises  yield  strength  Y  (or  y(p)),  the  calculated  zones 
of  crushed  concrete  were  inside  the  domain  of  the  plastic  zone,  Zo^^p.  Numerical  experi¬ 
ments  showed  (ref  8)  that  once  the  concrete  is  crushed  and  the  porosity  cp  is  "squeezed 
out,"  the  changes  in  the  density  of  the  crushed  concrete  become  negligible.  Thus, 
assuming  that  the  coordinate  z*  is  within  this  region  of  completely  crushed  concrete  (with 
Pt»cor7sf),  the  term  RhOf  equation  19  can  be  integrated,  which  results  in  the  following 


p{zo) »  \fhVz\z)  -  iptVz^o)  +  P(z )  +  2  -  f 

2  2  - - - - « 

A5_  ' - 


dr 


dz 


(20) 


r=0 


Rs 


In  general,  determination  of  the  parameters  Vz(z=Zo,z‘),  p(z‘ ),  Szz(z=Zo,z‘ ),  and  s„,r(z) 
requires  numerical  solution  of  the  system  of  equation  1  to  3.  However,  in  the  case  of  an 
axially  symmetrical  plastic  flow,  a  closed  form  analytical  solution  for  the  deviatoric  stress  Szs 
can  be  found  on  the  center  line  r=0.  In  the  case  of  the  von  Mises  yield  criterion,  equation  8 
to  10,  the  stress  deviatoric  tensor  is  given  by  the  following  flow  rule 


Sij 


>[7 


(21) 


where  is  the  strain  rate  tensor.  The  strain  rate  parameter  /  is  given  by  the  following 
equation 
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(22) 


'=  +el,  +  si)+i(yL+yi+Ti) 

Let's  evaluate  and  /  on  the  center  line  r=0.  Since  the  considered  plastic  flow  field 
is  axially  symmetric,  v^,r=v„^=Vr,^=v^,z=v^,^  =0.  It  follows  then  that  at  r=0 


dv. 
dr 


8«.  =  -^|  =0 


r=0 


Yoz  =0 
=  0 


(23a) 

(23b) 

(23c) 


Y.  =0  (23d) 

Substituting  equations  23  b  to  d  into  equations  21  and  22  results  in  the  following 

~  ^zr  ~  ®  (24) 


/=  \{K  +  +  ^l) 


(25) 


Assumption  that  the  target  material  is  incompressible  permits  significant  simplification 
of  the  equation  for  the  parameter  /.  Consider  a  perfectly  incompressible  material  for  which 
the  following  equation  holds 


+  See  +  e^=0 


(26) 


Substituting  equation  23a  into  equation  16  results  in: 


lr=0 


Thus, 


4 


St?  “• 


V3 


r=0 


(27) 

(28) 
(29) 
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It  is  interesting  to  note  that  equation  29  can  be  useful  for  numerical  validation  of  com¬ 
puter  codes  for  axially  symmetric  plastic  flow  analysis.  The  profiles  of  deviatoric  stresses 
Szz  shown  in  figure  6a  are  based  on  a  numerical  solution  of  equations  1  to  3  using  the  CALE 
code  (ref  1 )  and  are  in  excellent  agreement  with  the  closed  form  analytical  predictions, 
equation  29. 

In  the  case  of  the  constant  yield-strength  model,  equation  29  results  in 


ASz2=  0 


(30) 


In  the  case  of  the  pressure  dependent  yield-strength  model,  k{p)=Y{p)/\^ 

As^z  =  I  [ V(p(zo))  -y(p(z ))]  (31 ) 

Assuming  that  above  the  compaction  pressure  limit  pc,  the  strain-hardening  of  concrete  is 
insignificant  [i.e.,  dY/dp  is  relatively  small  (fig.  4)],  equation  31  results  in  A^2«0. 

Figure  7  compares  various  terms  of  equation  20.  Since  the  contribution  of  the 
deviatoric  stress  term  ASa  is  either  zero  or  negligible,  the  term  ASa  is  not  shown  in  the 
figure.  Unfortunately,  a  closed  form  analytic  solution  for  the  radial  derivative  of  the  shear 
stress  deviator  Sf2,r|r=o  is  not  available.  The  shear  stress  deviator  term  can  be  evaluated 
numerically  based  on  the  results  of  the  numerical  experiments.  The  calculated  profiles  of 
the  radial  derivative  of  the  shear  stress  deviator  s^.r  along  the  center  line  r=0  are  shown  in 
figure  6b.  Although  one  can  intuitively  surmise  that  changes  in  the  constitutive  model 
significantly  affect  the  deviatoric  stress,  the  results  of  the  analysis  shown  in  figure  7  indi¬ 
cate  that  the  relative  changes  in  the  magnitude  of  the  shearing  deviatoric  stress  term  Rs  are 
rather  negligible. 

Comparison  between  the  relative  magnitudes  of  terms  of  equation  20  shows  that  the 
terms  that  contribute  the  most  to  the  target/penetrator  interface  pressure  are  the  pressure 
term  p(z’)  and  the  "hydrodynamic"  term,  0.5pt[v2^(z‘)-Vz^(zo)].  As  discussed  in  the  previous 
section,  during  the  high  velocity  penetration,  the  pressure  in  the  target  material  fluctuates 
rapidly  and  within  rather  wide  limits.  These  rapid  and  large  fluctuations  prohibit  establish¬ 
ing  a  direct  relationship  between  the  target  constitutive  modeling  and  the  pressure  terms 
p(zo)  and  p(z‘),  unless  the  instantaneous  values  of  p(zo)  and  p(z‘)  (fig.  7)  are  time  averaged. 
In  contrast  to  the  pressure  term  p(z‘),  the  "hydrodynamic"  term  (which  represents  the 
kinetic  energy  of  the  target  material),  is  not  susceptible  neither  to  large  nor  to  rapid 
fluctuations,  which  makes  it  a  convenient  parameter  for  examining  effects  of  the  target 
constitutive  modeling. 
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In  order  to  establish  the  physical  meaning  of  the  "hydrodynamic"  term  Rh,  it  is  instruc¬ 
tive  to  compare  equations  19  and  20  with  the  equation  for  the  pressure  at  the  target/pene- 
trator  interface  given  by  the  modified  hydrodynamic  theory  of  penetration,  equation  1 3. 
Comparison  between  equations  13  and  19  indicates  that  these  equations  are  identical 
provided  that 


]_ 

2 


(32) 


[p(^  )  Szz(Zo)  -  SzziZ  )  - 


(33) 


The  definition  of  pressure  employed  in  equation  20  implies  that  the  pressure  is  nega¬ 
tive.  However,  according  to  equation  13,  p20,  which  makes  it  necessary  to  introduce  a 
minus  sign  on  the  right  hand  sides  of  equations  32  and  33.  It  follows  then  that  the  idealized 
penetration  rate  u,  per  modified  hydrodynamic  theory  of  penetration,  is  given  by 


1 


(34) 


Assuming  that  the  target  material  is  incompressible,  the  idealized  penetration  rate  u  is 
given  simply  by 


t/ « [  ^  [vi{zo)-vi{z)]\  2  (35) 

P/0 

According  to  the  modified  hydrodynamic  theory  of  penetration,  increases  in  the 
strength  factor  Rt  decrease  the  penetration  rate  u.  Similarly,  one  can  expect  that  the  con¬ 
stitutive  modeling  changes,  resulting  in  increases  in  the  target  resistance  to  penetration, 
have  to  decrease  the  idealized  penetration  rate  u.  This  hypothesis  can  be  supported  by 
the  following  results.  The  numerical  experiments  showed  that  in  the  case  of  the  constant 
yield-strength  model,  the  idealized  penetration  rate  was  u -0.1 23  cm/^is.  However,  after 
the  employed  constitutive  model  was  changed  to  the  pressure-dependent  yield-strength 
model  (which  increased  the  target  resistance  to  penetration),  the  idealized  penetration  rate 
was  decreased  to  u"=0.121  cm/ps.  Because  the  relative  change  in  the  magnitude  of  the 
penetration  rate  is  very  small,  the  accuracy  of  these  predictions  will  be  addressed  in  the 
following  section. 

The  calculated  trajectories  of  the  front  and  the  rear  of  the  penetrator  for  the  constant 
yield-strength  and  the  pressure  dependent  yield-strength  models  are  shown  in  figure  8. 
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Numerical  experiments  with  the  constant  yield-strength  model  showed  that  the  analysis  was 
not  affected  by  the  value  of  the  von  Mises  yield  parameter  Y,  which  resulted  in  identical 
trajectories  for  all  cases  studied.  However,  changing  the  employed  constitutive  model  to 
the  pressure-dependent  yield-strength  model,  resulted  in  a  slight  decrease  of  the  slope  of 
the  trajectory  of  the  penetrator  front,  while  the  trajectory  of  the  penetrator  rear  remained 
unchanged.  The  significance  of  this  result  is  as  follows. 

The  slopes  of  the  trajectories  of  the  penetrator  front  and  the  rear  represent  the 
penetration  rate  u  and  the  penetrator  rigid  body  velocity  v,  respectively.  According  to  the 
modified  hydrodynamic  theory  of  penetration,  the  penetrator  erosion  rate  is  given  by 

^  (36) 

where  /  is  the  instantaneous  length  of  the  penetrator.  Numerical  experiments  showed  that, 
after  the  concrete  constitutive  model  was  changed  to  the  pressure-dependent  yield- 
strength  model,  the  penetration  rate  u  was  decreased.  Since  the  rigid  body  velocity  v 
remained  unchanged,  the  absolute  value  of  the  penetrator  erosion  rate  \dlld1\  was 
increased  proportionally.  According  to  the  hydrodynamic  theory  of  penetration,  the 
penetration  is  accomplished  at  the  expense  of  the  penetrator  material.  Therefore,  the 
faster  the  penetrator  erodes,  the  faster  it  consumes  the  available  kinetic  energy  which, 
finally,  results  in  a  decrease  in  the  penetration.  This  is  precisely  the  reason  why  the 
pressure-dependent  yield-strength  model  resulted  in  a  decrease  in  the  penetration  depth. 
Thus,  it  is  reasonable  to  relate  the  notion  of  target  resistance  to  penetration  rate  u,  and  not 
to  the  pressure  at  the  target/penetrator  interface  p(zo),  nor  to  strength  related  parameters 
such  as  RsOr  ASzz. 


THE  EFFECT  OF  THE  TARGET  CONSTITUTIVE  MODELING 

Let  us  first  examine  the  effect  of  the  target  constitutive  modeling  on  the  axial  com¬ 
ponent  of  the  velocity  v^.  The  calculated  profiles  of  the  axial  component  of  velocity  v*  along 
the  center  line  r=0,  at  time  /=35  ps  after  the  impact,  are  shown  in  figure  9.  Note,  that  all 
curves  are  approximately  "parallel"  to  each  other,  and  that  increases  in  the  "spacing" 
between  the  curves  correlate  very  well  with  increases  in  the  target  resistance  to  penetra¬ 
tion.  Comparison  among  the  curves  shows  that  the  target/penetrator  interface  velocities 
Vz(zo)  differ  from  each  other  very  little.  Because  the  accuracy  of  assessing  values  of  Vz{zq) 
is  limited  by  the  resolution  of  the  computational  mesh,  the  effect  of  the  target  constitutive 
modeling  on  the  target/penetrator  interface  velocity  can  be  examined  more  conveniently  by 
employing  a  concept  of  the  idealized  penetration  rate  u,  given  by  equations  34  and  35. 
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As  discussed  in  the  preceding  section,  after  the  target  constitutive  model  was 
changed  from  the  constant  yield  strength  to  the  pressure-dependent  yield  strength,  the 
resulting  idealized  penetration  rate  was  reduced  by  a  relatively  small  value  of 
Ai/=i/-i/"=0.002  cm/ps.  The  accuracy  of  this  prediction  in  the  reduction  of  the  idealized 
penetration  rate  can  be  proved  by  results  of  numerical  experiments.  Assuming  that  the 
penetration  process  is  a  steady-state  one,  after  ^=35  ps,  the  difference  in  penetration  rates 
of  Au=0.002  cm/ps  produces  a  difference  between  the  axial  positions  of  the  target/pene- 
trator  interfaces  of  Azo«Au /=0.07  cm,  which  agrees  with  the  results  shown  in  figure  9. 

Note,  that  in  the  case  of  the  constant  yield-strength  model,  the  "spacing"  between  the 
curves  increases  consistently  with  increases  in  the  von  Mises  yield-strength  parameter  V, 
resulting  in  relatively  small  axial  "shifts"  of  the  Vz-curves.  In  the  case  of  the  pressure- 
dependent  yield-strength  model  the  "shift"  of  the  curve  is  the  largest,  which  implies  a 
relatively  significant  decrease  in  the  penetration  rate. 

One  should  also  note  that,  although  all  curves  shown  in  figure  9  are  approximately 
"parallel",  on  the  interval  [zo,z]  their  slopes  are  slightly  different,  which  is  related  to  the 
differences  in  the  target  resistance  to  penetration.  Compared  with  the  slopes  resulting  from 
calculations  with  the  constant  yield-strength  model,  the  change  to  the  pressure-dependent 
yield-strength  model  produced  a  rather  significant  decrease  of  the  "average"  slope  dvjdz 
(on  the  interval  [zo,z‘]).  Since  the  numerical  experiments  showed  that  regardless  of  the 
employed  constitutive  model  pt{z)  profiles  are  identical,  the  results  shown  in  figure  9  are  in 
complete  agreement  with  equation  34,  which  predicts  that  the  idealized  penetration  rate  u 
is  proportional  to  the  area  under  the  curve  /(z)=ptVz. 

in  the  proceeding  section  we  established  that  the  target  resistance  to  penetration  is 
directly  related  to  the  penetration  rate  u.  It  would  be  instructive  to  examine  the  structure  of 
the  term  Ru.  Rearranging  terms  of  equation  19  results  in  the  following  equation 


Rh  =  ^ptoU^=p(z^)-piz*)  +  5^(z*)-5^(Zo)  +  2[ 

2  ' - , - '  ' - , - ' 

Ap  -As^  ' — 


Since  changes  in  the  target  constitutive  modeling  have  a  negligible  effect  on  the  terms 
ASzzand  Rs,  it  follows  then  that  changes  in  the  penetration  rate  u  are  affected  mostly  by 
changes  in  the  pressure  gradient  Ap.  Thus,  the  target  resistance  to  penetration  should 
increase  with  decreases  in  the  pressure  gradient  Ap. 

As  discussed  in  the  preceding  section,  rapid  and  large  pressure  fluctuations,  which 
are  usually  present  in  the  numerical  analysis  of  the  high  velocity  penetration,  are  due  to  the 
interference  of  the  stress  waves.  Since  the  penetration  process  is  a  steady  state  one,  and 
since  for  the  conditions  considered  the  stress  wave  propagation  speeds  are  approximately 
constant,  the  spatial  gradients  of  the  pressure  waves  should  be  relatively  stable.  This  can 
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be  clarified  by  an  example  of  surface  waves  in  a  liquid.  For  a  stationary  observer,  fluctua¬ 
tions  in  the  height  of  surface  waves  are  rather  significant,  while  in  the  reference  frame  that 
moves  at  the  speed  of  the  wave,  the  shape  of  the  waves  stays  approximately  the  same.  In 
the  case  of  high  velocity  penetration,  numerical  experiments  show  that,  although  pressure 
profiles  may  slightly  vary  with  time,  their  shapes  are  similar  and,  therefore,  "average" 
spatial  gradients  dpidz  are  relatively  stable  regardless  of  the  employed  constitutive  model. 
Figure  10  shows  calculated  profiles  of  the  pressure  p  along  the  center  line  r=0,  at  time  f=35 
ps  after  impact.  Similar  to  the  decreases  in  slopes  of  the  Vz-curves,  the  shapes  of  pressure 
profiles  show  that  decreases  in  the  slopes  dpidz  correlate  very  well  with  increases  in  the 
target  resistance  to  penetration.  Thus,  the  results  of  the  numerical  experiments  are  in 
complete  agreement  with  the  predictions  of  equation  37,  which  states  that  the  penetration 
rate  is  proportional  to  the  pressure  gradient  Ap. 

This  general  tendency  (i.e.,  that  the  gradients  of  plastic  flow  field  variables  consis¬ 
tently  diminish  with  increases  in  the  target  resistance  to  penetration)  is  related  to  increases 
in  the  size  of  the  plastic  zone  discussed  in  the  preceding  section  and  shown  in  figure  2. 

The  nature  of  the  relationship  between  the  gradients  of  the  flow  field  variables  and  the 
extent  of  the  plastic  zone  is  as  follows:  numerical  experiments  show  that  the  changes  in 
the  target  constitutive  modeling  do  not  significantly  affect  plastic  flow  field  variables  neither 
at  the  target/penetrator  interface  Zo,  nor  at  the  elastic-plastic  boundary  Zep.  It  follows  then 
that,  since  inside  the  plastic  zone  Zo^^p flow  field  variables  change  monotonically,  once 
the  extent  of  the  plastic  zone  A=Zep-Zo  is  increased,  the  gradients  of  flow  field  variables 
dg{z)ldz-Ag{z)/A  are  decreased  proportionally  to  increases  in  A. 

A  brief  comment  is  needed  on  the  pressure  discontinuity  in  the  vicinity  of  z=z*  (fig.  10). 
As  discussed  in  the  previous  section,  the  employed  numerical  model  attempted  to  model 
the  behavior  of  the  break  gages.  The  results  of  numerical  experiments  show  that  during  a 
high  velocity  penetration,  the  behavior  of  thin  copper  sheets  surrounded  by  "crushable" 
concrete  is  similar  to  the  behavior  of  a  flexible  membrane  restraining  a  compressible 
continuum.  In  agreement  with  this  analogy,  the  numerical  analysis  shows  that  across  the 
thickness  of  a  break  gage  the  pressure  changes  very  rapidly,  which  results  in  a  "jump"  from 
the  pressure  build-up  "inside"  the  zone  of  completely  crushed  concrete  to  the  "ambient" 
pressure  "outside"  of  this  zone.  Small  differences  in  the  axial  positions  of  the  break  gages 
are  due  to  the  differences  in  the  target  resistance  to  penetration.  The  higher  the  target 
resistance  to  penetration,  the  smaller  was  the  displacement  of  the  break  gage. 

At  this  point,  it  is  appropriate  to  address  experimental  results  on  measuring  the  size  of 
the  plastic  zone  in  the  front  of  the  penetrator.  The  technique  for  measuring  the  plastic  zone 
size  in  the  front  of  the  penetrator  is  based  on  comparison  between  the  calculated  trajectory 
of  the  penetrator  front  and  the  experimental  data  from  trigger  times  of  the  break  gages. 

The  results  of  this  comparison  are  shown  in  figure  8.  Assuming  that  "breaks"  in  gages 
occurred  approximately  at  the  time  when  the  elastic-plastic  boundary  reached  the  gages, 
comparison  between  the  linear  regression  of  the  experimental  data  and  the  trajectory  of  the 
penetrator  front  results  in  the  size  of  the  plastic  zone  of  Aoxp«1 .3  cm.  This  value  of  the 
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plastic  zone  size  agrees  with  the  predictions  based  on  calculations  with  the  constant  yield- 
strength  model  [A«1 .3  cm  (fig.  2a)]  but  disagrees  with  the  predictions  resulting  from  the 
analysis  with  the  pressure-dependent  yield-strength  model  [A«1 .7  cm  (fig.2b)]. 

As  shown  in  figure  8,  the  trigger  time  data  are  scattered,  which  may  be  due  to  a 
number  of  reasons  including  inhomogeneities  in  concrete,  the  uncertainty  in  the  mode  of 
the  gage  failure,  etc.  It  should  be  noted  that  the  deviation  of  the  trigger  time  data  (with 
respect  to  their  linear  regression)  is  approximately  of  the  same  order  as  the  discrepancy 
between  Aexp=1 .3  cm,  A=1 .3  cm,  and  A=1 .7  cm.  Therefore,  with  the  present  degree  of 
accuracy,  this  technique  can  not  be  applied  for  proving  the  validity  of  either  of  the  consti¬ 
tutive  models  studied,  based  on  the  experimental  data  of  the  plastic  zone  size.  Previously, 
the  break  gage  technique  was  successfully  used  for  validation  of  the  concrete  equation  of 
state  modeling,  and  showed  that  the  "porous"  equation  of  state  model  is  the  most  realistic 
for  representing  the  elastic-plastic  flow  in  concrete  medium  resulting  from  the  high  velocity 
penetration  (ref  8).  Extension  of  this  technique  to  problems  related  to  the  concrete 
constitutive  modeling  would  require  a  number  of  improvements  which  ought  to  reduce  the 
amount  of  scatter  in  the  trigger  time  data.  This  can  be  accomplished  in  a  number  of  ways 
such  as  improving  the  homogeneity  of  concrete  (in  the  vicinity  of  the  break  gages),  refining 
the  gage  design  (to  exclude  the  uncertainty  in  the  gage  failure  mechanism),  etc. 

Another  important  penetration  parameter  that  is  significantly  affected  by  the  target 
constitutive  modeling  is  the  radial  component  of  velocity,  w.  Profiles  of  Vrat  a  distance  of 
approximately  0.7  cm  in  the  front  of  the  target/penetrator  interface,  at  ^=35  ps  after  impact, 
are  shown  in  figure  1 1 .  Comparison  between  these  curves  shows  that  with  increases  in  the 
target  resistance  to  penetration,  the  radial  component  of  the  velocity  is  consistently 
decreasing.  The  effect  of  this  parameter  on  the  formation  of  the  crater  is  as  follows. 

According  to  the  cavity-expansion  theory  (e.g.,  ref  26),  the  motion  of  the  target 
material  can  be  viewed  as  a  continuous  expansion  of  a  series  of  imaginary  cavities  which 
initially  have  zero  radius  and  then  expand  to  a  size  that  permits  passage  of  the  penetrator. 
When  the  impact  velocities  are  sufficiently  high,  the  penetrator  nose  deforms  into  a 
characteristic  mushroom-like  shape.  Behind  the  mushroom,  the  hole  continues  to  grow 
under  the  inertial  forces  that  are  due  to  the  radial  motion  of  the  target  material.  The  mag¬ 
nitude  of  these  forces  is  determined  by  the  intensity  of  the  target  material  flow  at  the  pene¬ 
trator  mushroom.  The  smaller  the  radial  component  of  the  velocity,  the  smaller  are  the 
inertial  forces,  and,  therefore,  the  smaller  are  the  resulting  crater  diameters. 

Another  important  parameter  that  determines  the  final  hole  diameter  is  the  size  of  the 
plastic  zone  behind  the  penetrator  mushroom.  Numerical  experiments  showed  (fig.  5)  that 
after  the  target  constitutive  model  was  changed  to  the  pressure-dependent  yield-strength 
model,  the  extent  of  the  plastic  zone  in  the  target  was  increased  not  only  in  front  of  the 
penetrator  but  also  behind  the  penetrator  mushroom.  Increases  in  the  plastic  zone  size 
increase  the  degree  of  the  viscous  dissipation  of  the  energy  (which  through  plastic  work  is 
converted  to  heat)  stored  in  the  expanding  hole.  Thus,  the  increase  in  the  plastic  zone  size 
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combined  with  the  decrease  in  the  radial  velocity  profile  resulted  in  a  significant  decrease 
in  the  final  crater  diameter  shown  in  figure  2. 

CONCLUSIONS 

1 .  Analysis  of  the  time  average  pressure  at  the  interface  between  the  penetrator 
and  the  target  shows  that,  in  an  agreement  with  the  hydrodynamic  theory  of 
penetration,  the  interface  pressure  increases  with  increases  in  the  "target 
strength,"  a  tendency  that  may  be  intuitively  associated  with  the  effect  of  the  von 
Mises  yield  strength  k  or  k{p).  However,  calculations  show  that  in  the  case  of 
the  pressure-dependent  yield  model  (the  "average"  yield  is  approximately  1 .6 
Kbar),  the  time  averaged  pressure  at  the  projectile/target  interface  is  higher  than 
the  pressure  resulting  from  the  1 .8  Kbar  constant  yield  model.  It  follows  then 
that  the  target  resistance  to  penetration  is  not  related  simply  either  to  the 
pressure  or  to  the  yield  strength  alone. 

2.  The  size  of  the  plastic  zone  in  front  of  the  penetrator  was  shown  to  be  directly 
related  to  target  resistance  to  penetration.  The  numerical  experiments  with 
different  constitutive  models  showed  that,  depending  on  the  particular  structure 
of  the  employed  model,  the  size  of  the  plastic  zone  A  was  either  significantly 
affected  by  the  model  or  the  effect  was  negligible.  In  particular,  all  numerical 
experiments  with  the  constant  yield-strength  model  resulted  in  approximately  the 
same  value  of  A«1.3  cm,  regardless  of  the  value  of  the  yield  strength  parameter 
/.  However,  a  change  in  the  target  constitutive  model  from  the  constant  yield- 
strength  model  to  the  pressure-dependent  resulted  in  a  very  significant  30% 
change  in  the  size  of  the  plastic  zone  (A=1 .7  cm). 

3.  Analysis  of  the  results  of  the  numerical  experiments  showed  that  the  penetration 
rate  u  is  the  principle  parameter  that  is  directly  related  to  the  target  resistance  to 
penetration.  Since  the  rigid  body  velocity  v  is  not  significantly  affected  by  the 
changes  in  u,  the  penetrator  erosion  rate  increases  proportionally  to  decreases 
in  the  penetration  rate.  The  faster  the  penetrator  erodes,  the  faster  it  consumes 
the  available  kinetic  energy,  which  ultimately  results  in  a  decrease  of  the  pene¬ 
tration. 

4.  Detailed  analysis  of  the  plastic  flow  field  at  the  center  line  showed  that  the 
penetration  rate  u  is  directly  proportional  to  the  pressure  gradient  Ap.  Numerical 
experiments  showed  that  changes  in  the  target  constitutive  modeling  do  not 
significantly  affect  plastic  flow  field  variables  neither  at  the  target/penetrator 
interface  Zo,  nor  at  the  elastic-plastic  boundary  Zep.  Given  that  inside  the  plastic 
zone  Zo^:^ep  flow  field  variables  change  monotonically,  once  the  extent  of  the 
plastic  zone  A=Zep-Zo  is  increased,  the  gradients  of  flow  field  variables 
dg{z)ldz-^Ag{z)IA  are  decreased  proportionally  to  increases  in  A.  Therefore,  the 
larger  the  size  of  the  plastic  zone  A,  the  slower  is  the  penetration  rate  u,  and  the 
higher  is  the  target  resistance  to  penetration. 
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)f  the  cx3mputational  domain;  spherical 
circular  cylinder  of  concrete  with  veloci 


□  Experimental 


Figure  2 

Resulting  hole  profiles  employing  (a)  the  constant  yield-strength  model  and  (b)  the 
pressure  dependent  yield-strength  model.  Experimental  data  referencel 
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Figure  3 

History  of  pressure  p=p(f)  at  the  target/penetrator  interface  at  r=0  for 

various  constitutive  models 
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Figure  4 

Empirical  yield  strength  function  Y-Y{p)  for  the  pressure  dependent 
yield-strength  model  (ref  16) 
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Figure  5 

Structure  of  the  elastic-plastic  flow  field  in  target  (time  t=35^,  y=0.062) 
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Penetrator/target  interface 
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Figure  6 

Calculated  profiles  of  (a)  Sa  and  (b)  Saj  along  the  center  line  r=0  for  the  constant  yield-strength  and 

the  pressure  dependent  yield-strength  models  (f=35  ps) 
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Figure  7 

Comparison  between  the  instantaneous  values  of  various  terms  comprising  the  pressure 

at  the  target/penetrator  interface  {t=35  ps) 


Figure  8 

Trajectories  of  the  front  and  the  rear  of  the  penetrator  for  the  constant  yield-strength  and  the 
pressure  dependent  yield-strength  models.  Experimental  data  reference  8. 


Figure  9 

Calculated  profiles  of  the  axial  component  of  velocity  Vj  at  the  center  line  r=0  (/=35  ps) 
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Figure  1 0 

Calculated  profiles  of  pressure  p  along  the  center  line  r=0  (f=35  ps) 


0.030 


stI/iuo  ‘m  Aiiooioa  |B|PBu 


32 


Figure  1 1 

Profiles  of  the  radical  velocity  component  Vr  versus  radial  coordinate  r(z-5.2  cm,  f=35  ps) 
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